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Abstract 

Using the basic concepts of chain by chain method we show that the 
symplectic analysis, which was claimed to be equivalent to the usual 
Dirac method, fails when second class constraints are present. We 
propose a modification in symplectic analysis that solves the problem. 
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1 Introduction 



There are some attempts to study a constrained system in the framework of 
first order Lagrangian [1, 2]. The coordinates appearing in a first order La- 
grangian are in fact the phase space coordinates. The Euler-Lagrange equa- 
tions of motion of a first order Lagrangian in an ordinary (non-constrained) 
system are the same as the canonical equations of motion. The kinetic term 
in a first order Lagrangian constitutes of a one-form whose exterior deriva- 
tive appears in the equations of motion. The resulted two-form, called the 
symplectic tensor, is singular for a constrained system. If the system is not 
constrained, usually the inverse of the symplectic tensor exists and provides 
the fundamental Poisson brackets (we exclude degenerate systems discussed 
in [17, 18] in which the symplectic tensor may have a lower rank in some 
regions of the phase space). 

The properties of a constrained system can be determined by trying to 
overcome the singularity of the symplectic tensor. Faddeev and Jackiw [3] 
used the Darboux theorem to separate canonical and non-canonical coor- 
dinates. They solved the equations of motion for non-canonical coordinates 
either to decrease the degrees of singularity of the symplectic tensor or to 
find the next level constraints. 

Then using a special system of coordinates, the authors of [4] showed 
that the Faddeev- Jackiw approach is essentially equivalent to the usual Dirac 
method [5]. In a parallel approach, known as symplectic analysis [6, 7, 8, 9] 
one extends the phase space to include the Lagrange multipliers. In this 
approach the consistency of constraints at each level adds some additional 
elements to the symplectic tensor. In other words, the kinetic part of the 
(first order) Lagrangian is responsible to impose the consistency conditions. 

The important point in most papers written in Faddeev- Jackiw method or 
symplectic analysis is that they often show their results for the constraints 
in the first level and then deduce that the same thing would be repeated 
at any level. However, the whole procedure of studying the singularities of 
symplectic tensor, demonstrates some global aspects. For example, some 
questions that may arise are as follows: 

What happens, after all, to the symplectic tensor? Is it ultimately singu- 
lar? How many degrees of singularity may it have? What is the relation of 
ultimate singularities with the gauge symmetries of the system? and so on. 
In [10] we showed that the symplectic analysis gives, at each step, the same 
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results as the traditional Dirac method (in the framework of level by level 
approach). The symplectic analysis may also be studied in the framework of 
chain by chain approach[ll] to obtain the Dirac constraints. 

Meanwhile, some recent observation [12] shows that in some examples the 
result of symplectic analysis and the well-established method of Dirac are 
not the same. This creates serious doubt about the validity of the symplectic 
analysis. Therefore, it is worth studying the origin of the difference between 
this approach and that of Dirac [5]. This is the aim of this paper. In the 
next section we first review the basic concept of symplectic approach as given 
in [10]. As we will show the symplectic analysis is equivalent to a special 
procedure in Dirac approach in which one uses the extended Hamiltonian at 
each level of consistency. In section(3) we will show that in the framework of 
Dirac method one is not allowed to use an extended Hamiltonian when there 
exist second class constraints. The important point to be emphasized is that 
this result can be understood more clearly in the framework of chain by chain 
method. In section(4) we show that for a one chain system with second class 
constraints the symplectic analysis as proposed in the literature fails. This 
result can be simply generalized to the general case of a multi-chain system. 
When recognizing the origin of the problem, we give our prescription to solve 
it in section(5). Finally in section(6) we give an example. 

The last point to be noticed is that the problem would not show itself 
for systems with two levels of constraints. As we will show, this is the case 
for second class systems with at least four levels of constraints. That is the 
reason for the fact that the problem does not appear if one considers just 
first level of constraints. 



2 Review of symplectic approach 

Consider a phase space with coordinates y l (i = 1, . . . , 2K) specified by the 
first order Lagrangian 

L = a i (y)y i -H(y) (1) 

where H(y) is the canonical Hamiltonian of the system. The equations of 
motion read 

fijy j = ihii (2) 
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where di = and the presymplectic tensor is defined as 

fij = didjbj) - djOi(y). (3) 

We denote it in matrix notation as /. This matrix is invertible for a regular 
system. Let f ij be the components of the inverse, From (2) we have 

y* = {y\H}, (4) 

where the Poisson bracket { , } is defined as 

{F(y),G(y)} = rd i Fd j G. (5) 

If / is singular, then using the Darboux theorem, as shown in [3], one can 
choose the independent coordinates (y' a , \ l ) such that 

L = a'j"* - \%(y') - H{y') (6) 

where f' aj3 = d a a'p — dpa' a is invertible. This shows that one can consider a 
system with a singular tensor fy, as a regular one described by 

L = a'J' a ~ H{y>) (7) 

together with by the primary constraints $i(y'). In other words, without 
losing the generality one can assume that one is at first given the first order 
Lagrangian (1) with a regular presymplectic two-form (3), and then the set 
of primary constraints — 1, • • • , M) are applied to the system. In this 

way the system is described by the Lagrangian 

L = a t y* - - H(y) (8) 

in the extended space {y\ A M ). The equations of motion (2) should be replaced 
in matrix form by 

' y\ ( dH \ 



f 












(9) 



which is equivalent to Eq. (2) together with the constraint equations $ ft = 
(/' 1.---.M). 

Now one should impose the consistency conditions <&ft = 0. To do this, 
one should extend the space to include new variables rf and add the term 
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j^t&W (or equivalently — to the Lagrangian (8). This leads in the 
extended space (y, A, 77) to the equations 



V 



/ 





A 











-A 





, 



A 



/ OH 



where the elements of the rectangular matrix A are given by 



A^i djQ^ ^ . 



(10) 



(11) 



However, nothing would be lost if one forgets about the variables \^ and 
reduces the system to the Lagrangian 



= a l f-i ] ^ -H(y). 
This leads to the symplectic two-form 



/ 


A 


-A 






(12) 



(13) 



in the (2K + M) dimensional space of variables Y = (y\r]^). It should be 
noted that the Lagrangian in Eq. (12) is the same as Eq. (8) in which A M 
is replaced by if . This means that the derivatives if have the same role as 
Lagrangian multipliers corresponding to primary constraints in the total 
Hamiltonian 

H T = H + \^£K (14) 

In other words, if some of ^'s are found by the dynamical equations of 
the system, then the corresponding Lagrange multipliers are obtained. In 
Dirac approach [14] this would be the case if there exist some second class 
constraints. 

The equations of motion due to the Lagrangian L* 1 ) can be written in 
matrix notation as 

FY = OH. (15) 
Using operations that keep the determinant invariant, it is easy to show that 



detF = det 



/ 







Af-'A / 



(det/) (det Af-'A). 



(16) 
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Since det / 7^ 0, F would be singular if C = A f 1 A is singular. Using (5) 
and (11) we have 

= (17) 

Suppose rank(C) = M" where M" < M. This means that F possesses 
M' = M — M" null-eigenvectors. One can, in principle, divide $^'s in two 



sets $iV and such that 



» C^v», det C^v» 7^ 0. 



(18) 



where the weak equality symbol ~ means equality on the surface of the 
constraints already known (here, the primary constraints). The matrix A 
can be decomposed to A' and A" such that 



i) 



A^ = d^>. 



(19) 



Accordingly the symplectic tensor F can be written as 







A" 


A' 




-A" 










\-A> 









(20) 



Consider the rectangular matrix 



{A'f-\0,l) 



(21) 



which has M' rows and IK + M columns. Using (18) one can show that its 
rows are left null-eigenvectors of F. Multiplying (21) with the equations of 
motion (15) gives the second level constraints as 



On the other hand, F in (20) has an invertible sub-block 



mv 



/ 


A" 








(22) 



(23) 
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with the inverse 



F~ l = 
mv 



j-i _ j-^A"C"- x N'f- x 


-f-^AI'C"- 1 




c"- 1 



(24) 



This can solve the equations of motion (15) for variables rf' to give 

" -C^"{^,H} (25) 



if 



where C^" u " is the inverse of C^i^,. Inserting this in the Lagrangian (12) 
gives 

= a t {y)tf - ^' ^ - H^{y) (26) 

where 

H^=H-{H,^}}C^"^. (27) 

In this way a number of Lagrange multipliers corresponding to the sec- 
ond class constraints are derived whose effect is only replacing the canonical 
Hamiltonian H with FfW. Now we can forget about them and suppose that 
we are given the primary constraints $^and the second level constraints <j)W . 
Next, we should consider the consistency of &W and add the term — r)% &W to 
the Lagrangian LW. Renaming the previous r/^'s as rji, the new Lagrangian 
would be 

£ (2) = a t (yW - f,^ - r£*f - H^(y) (28) 
this gives the symplectic two-form 





( f 


AW 


AW 


FW = 


-AW 










{ -AW 









(29) 



in the space (y, rjx, r] 2 ). Assuming that the composed matrix A = (A^°\ AW^ 

FW has the same from as (13). One should again proceed in the same way 
to find the null-eigenvectors as well as the invertible sub-block of FW , The 
process goes on in this and the subsequent steps as explained in more detail 
in [10]. 

The important point to be emphasized is that the Lagrangian 



L {n) = a % {yW-Y.€^ ] - H{n \y) 



(30) 



k=i 
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at the n-th level, say, is equivalent to a system with extended Hamiltonian 

n 

R (n) = Hi n-1) + £ ^k) (31) 
fe=l 

at that level. In other words, the symplectic analysis is equivalent to the 
Dirac approach in the context of level by level method provided that at 
each level one adds the new constraints with the corresponding Lagrange 
multipliers to the Hamiltonian. In fact this slight difference with the standard 
Dirac method may lead to some difficulties as we will see in the following 
section. 



3 The problem with extended Hamiltonian 

The extended Hamiltonian formalism is well-known in the context of first 
class constraints [13, 14]. In fact, it can be shown that the dynamical equation 

g = {g,H E }, (32) 

leads to the correct equation of motion provided that g is a gauge invariant 
quantity. In Eq. (32) the extended Hamiltonian H E is defined as 

H E = H + A m $ m (33) 

where $ m are only first class constraints (primary or secondary). For a first 
class system, the extended Hamiltonian can also be used step by step during 
the process of producing the constraints . In other words, when all of the 
constraints are first class, there is no difference whether one uses <£> = {<£>, H?} 
or $ = {$, H E }. 

Now we show that the extended Hamiltonian formalism in Dirac approach 
is not suitable when second class constraints are present. We show this point 
for a system with only one primary constraint, i. e. a one-chain system in 
the language of chain by chain method. We remember that for such a system 
level by level and chain by chain methods coincide. 

Consider a system with the canonical Hamiltonian H(y) and one primary 
constraint <t^\ The total Hamiltonian reads 

H T = H + X^ 1) . (34) 
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Suppose the consistency of $^ leads to $ ( ^ 2 - ) = H^. Then $( 3 ) emerges 

as |$^ 2 \ h}, and so on. The iterative process that produces the constraints 
is described by 

$("+!) = {$(»),#}. (35) 

The above procedure progresses unless 1& n \Ht} m or jcjjW, cj^ 1 ) j ^ 
at the last step N. In the former case the constraints in the chain are first 
class, i.e. commute with each other [11]; while in the latter all the constraints 
are second class which means that the matrix 



(36) 



is invertible. In this case the Lagrange multiplier A would finally be deter- 
mined as 

X = {$(JV) ) $(!)}• ^ 

Using the Jacobi identity, it is shown in [11] that the matrix C nrn in Eq. (36) 
has the following form 



/ 




c 



In other words 










Q2(N-1) 



C 1N 
(J2N 



\ 



o c^- 1 ) 2 

ATI qN2 



C (N-1)(N-1) C (N-1)N 
qN(N-I) qNN 



(38) 



(39) 



{$« 5 $(i)|~0 if i + j<N. 
Moreover using the Jacobi identity one can show from (35)that 

» - {$(9 » . . • » (_l)(f -1) {$(f ), $(f +D} ^ 0. 

(40) 

Remember that iV is the number of second class constraints and necessarily 
should be even. 

Now suppose that in order to define the dynamics of the system at some 
level n, one wishes to use the extended Hamiltonian 



(41) 



k=i 
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If n < y then from (38) the consistency of the constraint $ (n - ) gives 

$(n) = | $ (n) ) ~ |$(n) ) (43) 



which by (35), is the same as $( n+1 ). However at level y + 1 the consistency 
of $(f +1 ) ; using iJ^ 2+1 ^ gives 



$(# +D = {$(f +D, h) + A f $(f )} . (43) 



As is apparent from (40) the above equation solves the Lagrange multiplier 
\n. There is no justification to keep as the next constraint 

<3>(t+ 2 ). In order to knit the second class chain up to the last element & N \ 
one is just allowed to use the total Hamiltonian (34). In other words, the 
second half of the chain can be derived if only the primary constraint 
is present in the corresponding Hamiltonian. As explained in the previ- 
ous section, using the standard symplectic analysis is equivalent to using 
the extended Hamiltonian formalism described above. So one should expect 
some contradiction in symplectic analysis when second class constraints are 
present. In the next section we will show the essence of this contradiction 
for a one chain system and propose a method to resolve it. 



4 Second class one -chain in symplectic anal- 
ysis 



According to the algorithm given in section 2, given the canonical Hamilto- 
nian H(y) and the primary constraint at the first step of consistency 
one should consider the Lagrangian (see 12) 



=a i y i -r ) i® {1) -H(y). 



(44) 



The equations of motion can be written in matrix form as 




(45) 



It is easy to see that 



(i (1) rM 



) 



(46) 
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is the null-eigenvector of the matrix 



F 













(47) 



Implying u 1 on both sides of (45) and using (5) gives the new constraint 

$( 2 ) = (48) 
Adding the term — fj 2 ^W to the Lagrangian (to perform consistency) gives 



L W = a .yi _ ^$(1) _ ^$(2) _ H ( y y 
The equations of motion are 

dH \ 



(49) 



/ 


AW 


AW 


\ 


( * \ 


-AW 










-AW 








/ 


V V2 J 







(50) 



o J 



Assuming j^^ 1 ), $( 2 )| pa 0, one can find the new null eigenvector 

u 2 =(i< a >/-\0,l). 



(51) 



Multiplying u 2 by (50) gives the new constraint $® = |$^ 2 \ and so on. 

Suppose one wishes to proceed in this way to find the constraints of 
the chain discussed in the previous section, i.e. the second class chain 
<f>W^ ■ ■ ■ , cJiW with the algebra given in (38-40). Suppose the above pro- 
cedure has been proceeded up to the step y + 1 where the equations of 
motion are 

/ 



/ 


AW 




a& +1 ) \ 


( * \ 


( dH\ 


-AW 










m 

















-i(f+i) 







o / 




V o / 



(52) 



Clearly no more null-eigenvector can be find. In fact adding the column 
and row corresponding to the constraint $(t +1 ) has increased the rank of the 
matrix F by two. This means that the equations of motion can be solved to 
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find rf~*} and f]^~ +1 \ There is no way in the context of symplectic analysis 
to proceed further to find the remaining constraints &~ +2 \ ■ ■ ■ , $>( N ) of the 
chain. This is really the failure of traditional symplectic analysis. In fact this 
is the reason why the symplectic analysis has failed in the example given in 
[12] (Particle in hyper sphere). We will discuss this example in section(6). 

What we showed here is the failure of symplectic analysis for a second 
class system with only one primary constraint (i.e. a one chain system). 
However, one can easily observe that for an arbitrary system with several 
primary constraints again the symplectic analysis would fail. The reason is 
that for such a system some of the constraints driven at level n, i.e. may 
have non vanishing Poisson brackets with constraints of previous levels while 
commuting with primary constraints. As we know from Dirac approach, in 
such a case the Poisson brackets of these constraints with Hamiltonian give 
the next level constraints. Meanwhile, a little care on symplectic analysis 
shows that in this number of Lagrange Multipliers corresponding to 

non-primary constraints would be determined and there is no way to find the 
next level constraints. In this way, we conclude that the symplectic analysis 
would fail whenever second class constraints emerge at third level or higher. 



5 How to solve the problem 

In this section we try to find a way to maintain the symplectic analysis by im- 
posing some modifications. The origin of the problem is the fact that $(t +1 ) 
has non- vanishing Poisson bracket with <3> ( -~- ) . As a result, the symplectic 
two- form on the left hand side of Eq. (52), i.e. 



/ 


AW 






-AW 

























o J 



(53) 



does not possess a new null-eigenvector. If one could consider the vector 

u<T+i) = (^+ 1 )/- 1 ,0 J ,... J 0,l) > (54) 
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as a null-eigenvector, then by multiplying -u^ +1 ) on the right hand side of 
(52), one would obtain the next constraint as 



$(f+ 2 ) = {$(f+ 1 ), J ff} 



(55) 



To reach this goal one should truncate those columns of F which are 
located after AW . In other words, instead of F in Eq. (53) one should 
consider the rectangular matrix 



( f 



F = 



- A w 










(56) 



J 



Clearly «(t +1 ) in Eq. (54) is the null-eigenvector of F. It is obvious that 
if one does the same thing in the subsequent steps, one can produce all the 
remaining constraints of the chain, i.e. • • • , ($>( N \ In the last step the 

chain terminates, since 

But what is the justification to find the null-eigenvectors of F, i.e. the 
truncated F. In fact using Eq. (5) the set of equations 



( f 



-AW 



\ -AW 



AW \ 



is equivalent to 



Vi = { yi ,H + rii$W} 
$0) = o 





i = 
3 = 



(57) 



2K 
, N. 



(58) 



Remembering that fji has the same role as the Lagrange multiplier Ai corre- 
sponding to the primary constraint we see that Eq. (58) is the correct 
equation of motion 

y i = {y l ,H T }. (59) 

On the other hand, it is easy to sec that the equations of motion resulting 
from Eq. (52) can be written as 



Vi = {Vi, h e} 



(60) 
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where He contains all derived constraint (including second class ones). In 
fact as we explained before, the correct equations of motion are (58) and not 
(60). 

Therefore, if one wishes to proceed in the context of symplectic analysis, 
one should consider Eq.(57) instead of Eq.(52) 



6 Example 

Consider the Lagrangian 



(61) 



where q = (gi, 
Hamiltonian is 



q n ). The primary constraint is P v . The corresponding 
H = ip 2 - v (q 2 - l) (62) 



where p = (pi, • • • ,p n )- I n the usual Dirac approach, using the total Hamil- 
tonian H T — H + XP V , the consistency of = P v gives the following chain 
of constraints 

= P v 
$(2) = q 2 _ l 

$( 3 ) = 2q.p 
$( 4 ) = 2 (p 2 + 2vq 2 ) 



(63) 



As is apparent, $ < - 4 - ) and $ < - 3 - ) are conjugate to $ < - 1 - ) and $ ( ^ 2 - ) respectively. It is 
worth remembering that although $^ 3 ^ is second class, when reaching at third 
level, the process of consistency should not stop, i. e. it should be proceeded 
one level more to find $ ( - 4 ^ which is conjugate to the primary constraint 
In the symplectic approach the corresponding first order Lagrangian is 



L = pq + P v v - ^p 2 + v (q 2 - l) - XP V . 
This gives the singular presymplectic tensor 

F = 



f 












(64) 



(65) 
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where / is the usual (2n + 2) x (2n + 2) symplectic tensor: 



/ 














(66) 



The equations of motion for y % = (q, v, p, P v , A) are fayi = diH T where 
H T = H + XP V . Clearly this gives the canonical equation of motion with 
Hamiltonian Ht, together with the constraint equation P v = 0. Adding the 
consistency term —f\\P v to the Lagrangian (64), where r\\ is a new variable 
and forgetting about the term proportional to A (which just reproduces the 
primary constraint) one finds 



L« = pq + P v v - rjiPv ~ \v 2 + v (q 2 - l) . 
This gives the equations of motion 

diH 

where Y l = (q, v, p, P v , 771) . In the matrix form we have 



(67) 



(68) 



(69) 



where = (0,0,0,1). Here, bold zero(0) means a row vector with n 
zero components. Clearly = (0, — 1, 0, 0, 1) is the left null-eigenvector of 
F^. Multiplying the equations of motion (68) from the left by gives the 
constraint $^ 2 ^ = q 2 — 1. 





AM ' 


_i(i) 






In the next level we have the Lagrangian 

L( 2 ) = L-77 1 P,-7) 2 (q 2 -l) 



(70) 



written in the space Y % = (q, v, p, P v , r/i, r] 2 ). The corresponding symplectic 
tensor reads 



f 


AM 


AM 


-AM 








-AM 





, 



(71) 



where AM = (2q, 0, 0, 0). Clearly uM = (0, 0, 2q, 0, 0, 1) is the null-eigenvector 
of F 2 . Multiplying the equations of motion F^Y^ = diHr from the left by 
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u 



(2) 



gives the next level constraint &W = 2q.p. Again considering another 



variable 773, the third level Lagrangian would be 

= L- 7nP v - 7) 2 (q 2 - 1) - ?)3(2q.p) 
This gives the following symplectic tensor 



(72) 



( f 



F W 



-AW 

_i(2) 
_i(3) 



AO) AW AW \ 





















(73) 



where AW = (2p,0,2q, 0). Now the crucial point appears. That is, F^ 3 ) 
has no new null-eigenvector. In fact one expects that multiplying uW = 
(— 2q, 0, 2p, 0, 0, 0, 1) by the equations of motion due to LW gives the next 
constraint $( 4 ) = 2(p 2 + 2vq 2 ). However, it can be easily checked that 
U (3) F (3) _t 0. Moreover, uW (with one additional zero as the last element) is 
no more the null-eigenvector of F^ 3 ^. This means that adding the (2n + 5)th 
row and columns to FW has led to increasing the rank of FW by two. In other 
words, the equations of motion for f) 2 and f] 3 can be solved. Unfortunately 
without any modification there is no way to find the Lagrangian 

L< 4 > — L — r]^ - fj 2 (q 2 - l) - 77 3 (2q.p) - t) 4 (2 (p 2 + 2t;q 2 )) . (74) 
If we could find then we would be able to have 



F W = 



f 


AW 


AW 


AW 


AW 


-AW 














-Aw 














-Aw 














-Aw 











, 



(75) 



where AW = (8i>q, 4q 2 , 4p, 0). If we had somehow derived (74) and (75), 
then the singularity of symplectic tensor would completely disappear and 
771, • • -774 would be obtained. However, using the truncated symplectic tensor 
at the second step as 



/ 


AW 


-Aw 





-Aw 






(76) 
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and similarly F^ at the third level as 



^(3) 



V 



-AW 

_i(2) 



AW \ 



o J 



(77) 



makes it possible to introduce again and as the corresponding left 
null-eigenvectors of and F^ 3 \ respectively. This makes us able to find 
$( 4 ) as explained before. It should be noted that one can after all write the 
complete symplectic tensor F^ 4 ^. 

This example has also been discussed in [12], where some other reason is 
proposed as the origin of failure of the symplectic analysis. The same results 
as what we derived here can be found in every second class system possessing 
at least four levels of constraints. For example, one can study the simpler 
Lagrangian L = xy — z(x + y) as well as the more complicated example of 
bosonized Schwinger model in (1 + 1) dimensions [15, 16] given by 

C = l -d^dy + (<T - e^)d^A v - \f, u F^ + l -A^. (78) 

One can see that the main feature of the above calculations will more or less 
appear in all such examples. 
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